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We introduce a variational wave-function to study the polaron formation when the electronic 
transfer integral depends on the relative displacement between nearest-neighbor sites giving rise 
to a non-local electron-phonon coupling with optical phonon modes. We characterize the polaron 
crossover by analysing ground state properties such as the energy, the electron-lattice correlation 
function, the average phonon occupation and the quasiparticle spectral weight. Variational results 
are found in good agreement with numerical exact diagonalization of small clusters, and follow the 
correct perturbative result at weak coupling. We determine the polaronic phase diagram and we 
find that the tendency towards strong localization is hindered from the pathological sign change of 
the effective next-nearest-neighbor hopping. 



I. INTRODUCTION 

Significant electron-phonon (el-ph) interactions have been experimentally detected in many materials of wide inter- 
est, like manganites, 1 fullerenes, 2 , carbon nanotubes, 3,4 and cuprates. 5 In many of these cases, the el-ph interaction 
gives rise to polaronic features. A polaronic state results in fact when the electrons are strongly coupled to lattice 
distortion, therefore increasing their effective mass and leading to a state with low mobility. Increasing the el-ph cou- 
pling, the spatial extension of the lattice deformation decreases 6 and the polaron can vary its size from large to small. 
The single polaron problem of one electron interacting with the lattice degrees of freedom has been studied in detail, 
and allowed us to understand in the detail the physics leading to the formation of polaronic states. In particular it 
has been shown that the self-trapping process, which lead to the formation of polarons, is not a phase transition, but 
just a continuous crossover with no broken symmetry. 7 In the case of the Holstein model, 8 where quantum vibrations 
interact locally with the electrons, the crossover from large to small polaron has been extensively studied by several 
numerical techniques 9-15 and variational approaches. 16-18 In particular all the ground state properties of the Holstein 
model can be described with great accuracy by a variational approach 18 based on a linear superposition of Bloch 
states that describe weak and strong coupling polaron wave functions. 

The case of non-local interactions, that in general are also present in real materials, is much less understood. The 
coupling with acoustical phonons has been studied in order to explain the anomalous transport properties of non-local 
excitations, like solitons and polarons, in various ID systems. 19-22 In particular the tight- binding Su-Schrieffer-Heeger 
(SSH) model 19 was introduced to explain the transport properties of quasi one-dimensional polymers as polyacetylene 
where the CH monomers form chains of alternating double and single bonds. In this case the localization is due 
to a large shrink of two particular bonds and corresponding large hopping integral between the sites. As a result, 
the hopping between the two occupied sites and the surrounding ones is reduced resulting in a tendency towards 
localization. 

Our purpose here is to examine the single polaron formation in a model where both (Holstein) local and (SSH) 
non-local el-ph interactions are present. Due to the complexity of the model, we start the analysis using a perturbative 
approach that, although it can not capture the full multiphononic nature of the polaron, it has proved a remarkably 
useful tool in understanding the el-ph physics. 12,22 In particular we characterize the ground state properties of the 
system evaluating the energy, the electron-lattice correlation function and the quasiparticle spectral weight, in order 
to provide signatures of the polaron formation. In the limit when local el-ph interactions are much stronger or weaker 
than non-local el-ph interactions our model reduces to the standard Holstein model and to the SSH model with a 
dispersionlcss phonon spectrum, respectively. Then we start an accurate analysis of the non-local limit case, being 
this case not yet fully examined. Recently, the fully adiabatic regime of this model has been used to explain changes in 
carbon-nanotube length as a function of charge injection. 4 Furthermore the non-local case has been previously studied 
by one of us using exact diagonalization of small clusters up to four lattice sites, where an anomalous optical absorption 
has been identified. 12 In particular in Ref. [ 12] it is shown that the strong-coupling solution is characterized by an 



1 



unphysical sign change of the effective next-nearest-neighbor hopping which is missing when acoustical phonons are 
considered. 21 

In this work we improve the previous numerical analysis considering a six-site lattice and introduce a variational 
wave-function to investigate the thermodynamic limit of the system. The variational approach is based on a linear 
superposition of Bloch states that provide an excellent description of the lattice deformations on left and right bond 
of the polaron, respectively. The wave-function closely resembles a variational state previously proposed for the study 
of the Holstcin model 18 and for the SSH case it allows to describe polaron features in good agreement with exact 
numerical diagonalization results. The variational approach recovers the pathological behavior of the effective next- 
nearest-neighbor hopping pointing out that the non-physical region of parameters always prevents a strong localized 
solution. We also explicitly show that, when the phonon frequency is not really small, the considered non-local 
SSH interaction supplies a tendency to localize for the single carrier which can be more effective than the Holstein 
localization. 

The scheme of the paper is the following. In Sec. II we present the model and set down the notation. In Sec. Ill 
we discuss perturbative calculations showing the role of SSH cl-ph coupling with respect to the Holstcin contribution. 
Sec. IV is devoted to the presentation of the variational method in the limit of non-local el-ph interactions and its 
comparison with the exact diagonalization results. Sec. V reports our concluding remarks. 

II. THE MODEL 

In extremely general terms, the interaction between electron and harmonic lattice deformations is described by the 
Hamiltonian 

n = £ clt id ({ Xk })c j>rT + E|| + E + E *({**})4,<v. (!) 

where c ia (ci i<7 ) is the fermion creation (destruction) operator, a is the spin index, tij({xk}) is the electronic transfer 
integral for fixed lattice deformations {xk}, M is the ionic mass, Kij is the spring constant matrix, and ei({xk}) is 
the local energy of the electron. For small deviations from the equilibrium positions of the lattice we can approximate 
Uj({xk}) and ei({xk}) to be linear functions of the lattice displacements {xk} obtaining a general model with el-ph 
interactions. In particular, limiting the hopping to nearest-neighbor sites of a linear chain, we make the assumption 

t n +i, n ({xk}) = ~t + a(x n+ i - x n ) (2) 

typically employed for the derivation of the el-ph SSH interaction term, and 

e l {{x k }) = a\xi (3) 

generally used in order to deduce the local el-ph Holstein interaction. If spinless electrons and dispersionless Einstein 
phonons are considered, the model becomes 

H = -t^2(c]a +1 + c\ +1 Ci) + E aX i ai + HinU ( 4 ) 

i i 

where H int is 

Hint = gup y^(ctc» + i + c\ +1 c t )(a\ +l + a i+ i - a] - a t ) + giw ^ 4 c i( a i + a i)> ( 5 ) 

i i 

with a\ (en) the phonon creation (destruction) operator and wo the quantum of vibrational energy per site. The 
quantity g = a/ \J2MuJq is the SSH coupling that we mainly discuss in this work, while g± — ct\/ ^2Mwg is the 
Holstein local electron-phonon coupling. We study the coupling of a single electron to lattice deformations using units 
such that the lattice spacing a = 1 and h = 1. 

III. PERTURBATION THEORY 

Weak-coupling perturbation theory in the electron-phonon coupling has proved a remarkably useful tool in under- 
standing the el-ph physics. Besides the obvious ability to describe the weakly interacting regime, the perturbative 
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approach has in fact provided some guidelines to understand the conditions for polaron formation in the Holstcin 
model. More explicitly, the polaron crossover occurs around the coupling value for which the perturbative approach 
breaks down. 12 

Here we focus on the case of one electron in a one-dimensional chain. If the el-ph terms are smaller than both 
the hopping term and the bare phonon term (g, gi <C t, 1 and t = t/uj ), they can be treated as perturbations of the 
unperturbed Hamiltonian Hq = Hki n + H p h ■ 

The second-order correction to the energy of the ground state is given by 



AE(0) 



l + 2t- VT+ At 
t 2 



9l 



\J\ + At 



(6) 



while the perturbative correction to the free band Sk = — 2tcos(fc) is reported in Appendix A. 

We note that for fixed values of the coupling constants g and g\ the two contributions (SSH-like and Holstein) 
have different behaviors as functions of the inverse adiabatic ratio t. In particular, as shown in Fig. 1, the Holstcin 
contribution to the ground state energy is always lower (for g — gi) than the SSH one when t < t w with t w = 4 + 3\/2. 
In other words, when the phonon frequency are not really small, the SSH cl-ph coupling is more effective than the 
Holstein one. The reduced effect of the Holstein el-ph coupling when t is small (anti-adiabatic regime) pushes the 
polaron crossover to larger values of the coupling A = g 2 /2u t as the phonon frequency is increased. Actually, while 
A > 1 is the condition for the polaron crossover in the adiabatic regime f > 1, in the antiadibatic regime t <C 1, it 
has been shown that the crossover occurs when a 2 = g 2 /oJq ~ 1, i.e., for A ^> l. 12 > 14 > 15 Recently it has been shown 
that this important role of the degree of adiabaticity is not limited to the single polaron problem, but it also extends 
to finite dennsities. 23 

Since polaron formation is not a phase transition and occurs without symmetry breaking, different criteria can be 
established to define the crossover values of the coupling constants which mark the polaronic regime. In the following 
we compute some physical quantities which have been often introduced to characterize the polaron crossover. 

The average phonon occupation number N p h — jj(J2i a \ a i) lSi gi ven by 



N ph = T 



(1 + 2t) 
\J\ + At 



+ 2\ x t 



r (1 + 2?) 



(1 + At) 3 / 2 ' 



(7) 



where A = g 2 ui /2t and Ai = g 2 u /2t. From Eq. (7) it turns out that the phonon number, as the ground state energy, 
is more affected by the SSH coupling when luq exceeds a given value. In particular for t < 2 (i.e. Uo/t > 0.5) the SSH 
contribution is always higher than the Holstein one. 

Other quantities of great interest to characterize the polaron formation are the electron-lattice correlation functions. 
In particular we consider the correlation function Xi,s — ( c l c i( a l+s + a i+s)) between the electronic density on a site 
i and the lattice displacement on site i + 5, which measures the entanglement of lattice and electronic degrees of 
freedom typical of the polaronic state. After a Fourier transformation in the momentum space, at k — one has 



Xfc=o,<5=o 



25i 



Xfc=0,<5=l 



Xfc=0,<5=2 



-2<? 



1 f 

7 + ^|3 



At 



l + 2t 
\/l + At, 



2gi 



Vl + Ai 



+ 



l + 2t 
2t 2 



1 - 



l + 2t 
Vl + Ai 



(8) 



In Fig. 2 we plot the correlation function at nearest-neighbor (left) and next-nearest-neighbor (right) sites as functions 
of the inverse adiabatic ratio t, for fixed values of the couplings g = gi = 1. As expected, the value of the correlation 
function goes to zero for large values of t, but the behavior of the SSH-like contribution (dashed lines) is qualitatively 
different from that of the Holstein ones (dot-dashed lines). 

The last quantity we consider is the quasiparticle spectral weight Z(k) = (1 — ^q^^ \u>=e(k))~ 1 ; which measures 
the renormalization of the electron Green's function due to the el-ph interaction. The second-order perturbative 
self-energy E(fc,w) is given in Appendix A. Even if the polaronic regime cannot be attained within lowest-order 
perturbative approach, indications on the beginning of the polaronic crossover can be extracted from the spectral 
weight expression. In particular the polaron crossover is expected to be associated with a sharp reduction of this 
quantity as a function of the couplings. The expression of the inverse spectral weight at k — is given by 
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(1 + 2t) 



+ 2Aii 



- (l + 2t) 



(1 + At) 3 / 2 



(9) 



while the full momentum dependence of Z{k) is reported in Appendix A. As expected the spectral weight Z(0) is a 
monotonically increasing function of t, for a fixed value of the couplings. It is interesting to note that the reduction 
of Z(0) due to the SSH-like contribution is more relevant of the Holstein ones for t < 2, while in the adiabatic limit, 
i.e. for large value of t, it is very small and slow. 

Strictly speaking, perturbative calculations only correctly characterize the small coupling regime. In order to provide 
a better insight on the problem of the polaron formation in the systems with non-local interactions, in the following 
we focus on the SSH contribution and substantiate our results by analytic variational calculations and numerical exact 
data. 



IV. VARIATIONAL APPROACH VS. EXACT DIAGONALIZATION 



In this section we extend our analysis of the non-local SSH model to the whole range of el-ph couplings using two 
standard and well grounded techniques, a variational approach and exact diagonalization of small clusters. First we 
introduce the variational wave function. We consider translation-invariant Bloch states obtained by superposition 
of localized states centered on different lattice sites. 24 These wave-functions have been introduced in order to study 
the polaron formation within the Holstein model where they are able to fully capture the features of the Holstein 
polaron. 16,18 In this work we extend this kind of wave- functions to the SSH interaction model assuming 

\^ } ) = ^zZ eik ' n \^\n)), (10) 

* n 

where \^\n)) is defined as 

\^\n)) = e^ i) W + ^ ) ("- 1 ) + ^ ) (" +1 )]|0) p ^^(m)e ife -ct +ro |0) ei , (11) 



with the quantity U^\j) given by 



U<£\j) = -2= £[/g-(<zKe^ - h.c.]. (12) 



The phonon distribution function f^- (q) is chosen as 



(») 

1 + 2i/?£j [cos(k) - cos(fc + q)] 

with a^j and variational parameters. In Eq. (11), \0} p h an d |0) e ; denote the phonon and electron vacuum state, 
respectively, and the variational functions (f>^ (m) are assumed to be 

4>f{™) = £ ^\j)S m ,j, (14) 

j=-5 

where 7^(j) are variational parameters that take into account the broadening of the electron wave- function up to 
fifth neighbors. It is worth to note that traditional variational approaches to the Holstein polaron problem uses the 
localized state (11) where only the on-site operator uj^\n) is applied. Thus we introduce in the expression of the trial 
wave-function the nearest-neighbor displacement operators U^\n + 1) and uj?\n — 1), in order to take into account 
the dependence of the hopping integral on the relative distance between two adjacent ions. 

Reflecting the asymmetry of the SSH coupling (shrinking of the bond on which the electron is localized and stretching 
of the neighboring bonds), we also define two wave- functions that provide a very good description of the lattice 
deformations on left and right bonds of the polaron. Naturally the left and right directions are relative to the site 
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where the presence of the electron is more probable. Thus in Eq. (11) the apex i = L,R indicates the Left (L) and 
Right (R) polaron wave-function, respectively. The wave-functions L and R are related as follows 



/g(«) = -/S(?)<o 
/ffi-i(ff) = -/ffi-i(ff)>° 

/2 + i(<?) = -/S + i(9)>° 

<i>l R) (m)=<j>i L \-m). (15) 

All the variational parameters are determined by minimizing the expectation value of the Hamiltonian (4) with gi = 
on the states (11). Even though the wave- functions L and R describe correctly the lattice deformations of the left and 
right side of the polaron, respectively, the mean values of the Hamiltonian on these states are equal. So the relations 
(15) can be also viewed as those that leave unchanged the energy functional determined by one wave- function. 

These two wave-functions can be improved by increasing the extension of the phonon contributions in Eq. (11) and 
of the electron terms in Eq. (14). Furthermore, they are not orthogonal and the off-diagonal matrix elements of the 
Hamiltonian between these two states are not zero. This allows to determine the ground-state energy by considering 
as trial state the linear superposition 18 of the wave- functions R and L 

where \& k L ^} and are the normalized wave-functions L and R weighted by the coefficients A k and B k and 

Sk = H L) H R) ) (17) 

is the overlap factor. The wave-function (16) correctly describes the properties of the lattice deformations on both 
the sides of the polaron and we will find that it is in very good agreement with the results derived by the exact diag- 
onalizations on a chain of 6 sites. Furthermore the variational approach involves a number of variational parameters 
that does not depend on the length of chain, so it allows to study the thermodynamic limit of the system. 

The minimization procedure is performed in two steps. First the energies of the left and right wave-functions 
are separately minimized, then these wave-functions are used in the minimization procedure of the quantity E k = 
(ip k \H\tp k ) I '(V'felV'fe) with respect to A k and B k defined in (16). 18 Exploiting the equality 



(4 L) \HHi L) )^(4 R) \H\4 R) )=e k , (18) 



we obtain 



£fc — S k E kc — \E kc — S k e k \ 



^ = - —>cj->c ^ (19) 



where E kc = ($ k L ^ \H \& k R ^} is the off-diagonal matrix element, and \A k \ = \B k \. The matrix elements between the 

states an d tp k L ^ contained in Eq. (19) are reported in Appendix B. 

The total energy functional (19) is minimized with respect to the variational parameters and the optimal ground 
state energy is plotted in Fig. 3 for a six-site lattice and two different values of the inverse adiabatic parameter t. We 
also study the thermodynamic limit and find energy curves very close to those of the finite system. In order to test 
the validity of our variational approach (VA), we perform exact numerical calculations on small clusters by means 
of the Lanczos algorithm. We improve the previous exact diagonalization (ED) analysis of the model, investigating 
small clusters up to six sites. 12 As shown in Fig. 3, each variational and exact numerical curve exhibits a kink with 
increasing the el-ph coupling. We have checked that at these couplings the effective next-nearest-neighbor hopping 
changes sign opening an unphysical region of the parameters. The agreement between numerical data and variational 
approach is very good up to g values close to the unphysical transition. 

In order to characterize the polaron formation we also analyze the electron-lattice correlation function x%,5 defined 
in Sec. III. In particular in Fig. 4 we show the behavior of Xi,6 as a function of the SSH coupling for S = 0, 1,2 
and t — 2.5. As expected, variational results and exact numerical data always recover the perturbative values in 
the limit of small cl-ph coupling. Increasing g the monotonic behavior of the correlation function exhibits a kink, 
as the ground state energy. In particular the correlation function at next-nearest-neighbor (5 = 2) changes sign as 
the effective hopping, confirming the pathological behavior. At couplings where the ground state energy and the 
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correlation function show the kink, also the average phonon number is characterized by an anomalous behavior as 
shown in the bottom right panel of Fig. 4. 

In order to extract information on the values of g at which polaron crossover begins, before the opening of the 
unphysical region, we also investigate the behavior of the quasipartcle spectral weight Z(0). We find that increasing 
the el-ph coupling for fixed values of t, the spectral weight starts to drop but it never reaches a really small value 
before the unphysical sign change of the hopping occurs. Nevertheless we observe distinct signatures of the tendency 
towards localization, as shown in Fig. 5, where Z(0) is plotted as a function of g for the fixed value t — 2.5. 

We conclude our analysis collecting the obtained data in the phase diagram of Fig. 6. It is calculated from the 
position of the kink in the ground state energy obtained by means of the variational approach (diamonds) and the exact 
diagonalization (triangles). The agreement between the two methods becomes better moving towards the adiabatic 
limit. In analogy with the phase diagram obtained for the Holstein polaron, 18 we also mark a crossover region defined 
as the range of parameters for which Z(0) is less than 0.9. As shown in Fig. 6, we find that the considered SSH 
model does not present any marked mixing of electronic and phononic degrees of freedom, being the strongly coupled 
state prevented from the pathology of the model. As far as the fully adiabatic limit luq = is concerned, we verify 
that the crossover line joins onto the line for the transition to the unphysical region at the critical value A = 0.25, 
confirming the discussion in Ref. [ 12]. We finally notice that, as discussed in Ref. [ 12], both the crossover region 
boundary, and the instability line obtained by exact diagonalization are only weakly dependent on the adiabatic ratio, 
and that A is the relevant electron-phonon coupling regardless the value of t. This is a peculiarity of the SSH coupling 
with respect to the Holstein one, where the polaron crossover moves to large values of A as the phonon frequency 



In this work we discussed the features of one electron non-locally interacting with optical phonons in a discrete 
chain. We introduced a variational wave function to locate the crossover region for the transition between weak and 
strong localized polaron solutions. In particular we found that the pathological sign change of the effective next- 
nearest-neighbor hopping always precedes a stable strongly localized solution. Such an unphysical region of the model 
parameters does not occur in the case of acoustical phonons being the deformation linked to the particle extension 
along the entire chain. 21 However we have also shown that, for finite values of the adiabaticity parameter, when the 
phonon frequency is not really small, the non local (SSH) cl-ph interaction is more effective than the local (Holstein) 
one in reducing the mobility of the electron. Then our variational calculations are an interesting starting point to 
examine the complex problem of the polaron formation in a model where both local and non-local el-ph interactions 
are present. In particular we emphasize that the proposed variational wave function for the SSH limit can be slightly 
modified to be suitable for the treatment of the complex case where both interactions are present. Detailed future 
investigations in this direction are required. Finally we stress that the validity of our variational results is supported 
by an accurate analysis of exact diagonalization data on small clusters. The agreement between VA and ED data is 
good up to coupling values close to the unphysical region. 
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In the limit of small el-ph couplings, the perturbative second order correction AE(k) to the tight-binding free band 
energy is 
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Moreover, using the bare phonon and electronic Green propagators, the perturbative self-energy reads 
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APPENDIX B 



In this appendix we report the matrix elements between the states \ip^ > and \ip^ >. These quantities are 



involved in the calculation of the ground-state energy within the variational approach. We find 

(4 L) \4 R) )= E ^(-mi^N^N-^). 



where the phonon matrix clement Z^ L R \i — j) is defined as 
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The quantity e k = (4 k L) \H\4 k L) ) = (^H^IVO is easily derived using the matrix elements given above. 
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FIG. 1. Second order correction to the gound state energy, Eq. (6), as a function of the adiabatic inverse ratio t, for 
g = gi — 1 (solid line). The dashed line is the SSH-like contribution (g — 1 and g\ = 0), the dot-dashed line is the Holstein 
contribution (g — and g\ = 1). 
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FIG. 2. Left: Correlation functions at nearest- neighbor sites (on the left) and next-nearest-neighbor sites (on the right) at 
k — as functions of the inverse adiabatic ratio t, for g = <?i = 1 (solid line). The dashed lines show the SSH-like contribution 
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FIG. 3. Ground state energy E(0) as a function of the SSH el-ph coupling g for two different values of the inverse adiabatic 
ratio i = 2.5 (left) and i = 1 (right). Solid and dotted lines are obtained from the variational approach and the Lanczos data 
for a six-site lattice, respectively; perturbative curves (dot-dashed lines) are plotted for comparison. Symbols mark the kink 
values of the energy. 
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(top left), 5—1 (top right) and 5 = 2 (bottom left) as functions of the 
SSH coupling g for t = 2.5. Bottom right: Phonon number vs. g for the same value of t. Solid lines are obtained from the 
variational approach in the thermodynamical limit; dotted lines show Lanczos data; perturbative curves from Eqs. (8) and Eq. 
(7) with gi = (dot-dashed lines) are plotted for comparison. 
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FIG. 5. Spectral weight Z(0) as a function of the SSH el-ph coupling g for t = 2.5. The solid line is obtained from the 
variational approach in the thermodynamical limit; the dotted line shows Lanczos data; the perturbative curve from Eq. (9) 
with gi = (dot-dashed line) is plotted for comparison. 
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FIG. 6. Phase diagram for one electron in a six-site lattice. Triangles and diamonds correspond, respectively, to the couplings 
where the exact numerical ground state energy and the variational result have a kink. The dashed line indicates the boundary 
of the crossover region, where the spectral weight Z(0) is less than 0.9. 
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